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Recent experimental results on hadronic r decays into strange particles by the OPAL collaboration 
are employed to determine V us and m s from moments of the invariant mass distribution. Our 
results are V u3 = 0.2208 ± 0.0034 and m s (2GeV) = 81 ± 22MeV. The error on Vus is dominated 
by experiment, and should be improvable in the future. Nevertheless, already now our result is 
competitive to the standard extraction of V us from K e 3 decays, and it is compatible with unitarity. 



PACS numbers: 12.15.Ff, 14.60.Fg, 11.55.Hx, 12.38.Lg 
INTRODUCTION 

Already more than a decade ago it was realised that 
the hadronic decay of the r lepton could serve as an 
ideal system to study low-energy QCD under rather clean 
conditions [1]. In the following years, detailed investi- 
gations of the t hadronic width as well as invariant 
mass distributions have served to determine the QCD 
coupling a s to a precision competitive with the current 
world average [2, 3]. The experimental separation of the 
Cabibbo-allowed decays and Cabibbo-suppressed modes 
into strange particles opened a means to also determine 
the mass of the strange quark [4-12], one of the funda- 
mental QCD parameters within the Standard Model. 

These determinations suffer from large QCD correc- 
tions to the contributions of scalar and pseudoscalar cor- 
relation functions [1, 12-14] which are additionally am- 
plified by the particular weight functions which appear 
in the r sum rule. A natural remedy to circumvent this 
problem is to replace the QCD expressions of scalar and 
pseudoscalar correlators by corresponding phenomeno- 
logical hadronic parametrisations [4, 7, 9, 10, 15], which 
turn out to be more precise than their QCD counter- 
parts, since the by far dominant contribution stems from 
the well known kaon pole. 

Additional suppressed contributions to the pseu- 
doscalar correlators come from the pion pole as well as 
higher exited pseudoscalar states whose parameters have 
recently been estimated [16]. The remaining strangeness- 
changing scalar spectral function has been extracted from 
a study of S-wave Kir scattering [17, 18] in the framework 
of resonance chiral perturbation theory [19]. The result- 
ing scalar spectral function was also employed to directly 
determine m s from a purely scalar QCD sum rule [20] . 

Nevertheless, as was already realised in the first 
works on strange mass determinations from the Cabibbo- 
suppressed r decays, m s turns out to depend sensitively 



on the element V us of the quark-mixing (CKM) matrix. 
With the theoretical improvements in the r sum rule 
mentioned above, in fact V us represents one of the dom- 
inant uncertainties for m s . Thus it appears natural to 
actually determine V us with an input for m s as obtained 
from other sources [4]. 

Very recently, new results on the r branching fractions 
into strange particles have been presented by CLEO [21] 
and OPAL [22] . In addition, the OPAL collaboration also 
presented an update on the strange spectral function, 
previously known only from ALEPH [10]. Both, CLEO 
and OPAL found B[t~ — > A'~7r + 7r~// T ] to be significantly 
higher than the corresponding ALEPH result. The im- 
portant impact of these improved findings on the deter- 
mination of V U8 and m s will be investigated below. 



THEORETICAL FRAMEWORK 

The main quantity of interest for the following analysis 
is the hadronic decay rate of the r lepton, 



_ T[t — > hadrons v T (7)] 
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which experimentally can be decomposed into a compo- 
nent with net-strangeness R T ,s, and the non-strange part 
Rt.ns- Additional information can be inferred from the 
measured invariant mass distribution of the final state 
hadrons. The corresponding moments Rf , defined by [23] 
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can be calculated in complete analogy to R T = R®°. In 
the framework of the operator product expansion (OPE), 
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can be written as [1]: 



K l = 3S EW {(|K d | 2 + |K,| 2 ) (l + <5 fe ^) 



+ E(i^i 2 ^r+iK s i 2 € (D) )}- ( 3 ) 

The electroweak radiative correction Sew = 1-0201 ± 
0.0003 [24-26] has been pulled out explicitly, and <5 feZ (°) 
denotes the purely perturbative dimension-zero contri- 
bution. The symbols 5^°^ stand for higher dimensional 
corrections in the OPE from dimension D > 2 oper- 
ators which contain implicit suppression factors l/M^P 
[9, 12, 13]. 

The separate measurement of Cabibbo-allowed as well 
as Cabibbo-suppressed decay widths of the r lepton [10, 
21, 22] allows one to pin down the flavour SU(3)-breaking 
effects, dominantly induced by the strange quark mass. 
Defining the differences 
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many theoretical uncertainties drop out since these ob- 
servables vanish in the SU(3) limit. 



DETERMINATION OF Vus 

Employing the SU(3)-breaking difference (4), as a first 
step, we intend to determine V us . This approach re- 
quires a value for the strange mass from other sources 
as an input so that we are in a position to calculate 
5Rf from theory. In the following, we shall use the re- 
sult m s (2GeV) = 95 ± 20MeV, a value compatible with 
most recent determinations of m s from QCD sum rules 
[16, 20, 27] and lattice QCD [28-30]. The compilation of 
recent strange mass determinations is displayed in fig- 
ure 1. For comparison, in figure 1, we also display m s as 
obtained from our previous r sum rule analysis [4] for the 
ALEPH data, as well as this work analysing the OPAL 
data. 

Since the sensitivity of 8Rf to V us is strongest for the 
(0,0) moment, where also the theoretical uncertainties are 
smallest, this moment will be used for the determination 
of V us . Inserting the above strange mass value into the 
theoretical expression for SR T [4] , one finds 



SR Tth = 0.218 ±0.026, 



(5) 



where the uncertainty dominantly results from a vari- 
ation of m s within its errors. Employing the above re- 
sult in eq. (4), together with the experimental findings 
R t NS = 3.469 ± 0.014, R TiS = 0.1677 ± 0.0050 [22], as 
well as \V ud \= 0.9738 ± 0.0005 [31], we then obtain 

\V US \ = 0.2208 ±0.0033 cxp ± 0.0009th = 0.2208 ± 0.0034 . 

(6) 
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FIG. 1: Summary of recent QCD sum rule [4, 16, 20, 27] and 
lattice QCD [28-30] results for m s (2GeV). For ref. [30] the 
given error is only statistical. 



The first given error is the experimental uncertainty, 
dominantly due to R T ,s, whereas the second error stems 
from the theoretical quantity 5R Tt th- For the extraction 
of V us , even though the theoretical error on 5R T .th is 
12%, it represents only a small correction compared to 
Rt,ns/\V u( i\ 2 and thus its error is suppressed. The the- 
oretical uncertainty in <5i? Ti th will only start to matter 
once the experimental error on i? Tj s is much improved, 
possibly through analyses of the BABAR and BELLE t 
data samples. 

One further remark is in order. A sizeable fraction 
of the strange branching ratio is due to the decay 
t -> Ku T , for which OPAL used the PDG fit result 
B[t -> Kv T {-y)} = (0.686 ± 0.023)% [31]. However, this 
decay can be predicted employing its relation to the de- 
cay K — > /ii^(7), which theoretically is known rather 
well [32, 33]. Updating the numerics of refs. [32, 33], 
we then obtain B[t -> Kv T {-y)] = (0.715 ± 0.004)%, 
much more precise than the experimental value. Adding 
this result to the remaining strange branching fractions, 
one finds R TjS = 0.1694 ± 0.0049, which would lead to 
\V US \ = 0.2219 ±0.0034. 



STRANGE QUARK MASS 

Employing the above calculated value for V us , we are 
now in a position to determine the strange quark mass 
m s from the SU(3)-breaking difference of eq. (4). Exper- 
imentally, various (k, I) moments have been determined 
[22]. For low k, the higher-energy region of the experi- 
mental spectrum, which is less well known, plays a larger 
role and thus in this region the experimental uncertain- 
ties dominate the strange mass determination, whereas 
for higher k more emphasis is put on the lower-energy 
region, and there the theoretical uncertainties dominate. 
At present, the most reliable results for m s are obtained 
from the moments (2,0) to (4,0), and we shall only dis- 
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Parameter 


Value 


(2,0) 


(3,0) 


(4,0) 


m s (M T ) 




93.2 


86.3 


79.2 


Tykl 
U t,NS 


[22] 


+ 5.1 

-5.4 


+3.6 

-3.7 


+2.8 
-2.9 


H T,S 


[22] 


—30.9 

+23.3 


— 19.5 
+ 15.8 


— 13.9 

+ 11.6 


\v us \ 


0.2208 ± 0.0034 


+21.7 
-29.8 


+ 14.6 
-18.7 


+ 10.6 
-13.0 




z X LJ{ot B ) 
no O(af) 


—4.0 
+4.6 


— 5.3 
+6.5 


—6.1 

+7.8 


I 


1.5 
0.75 


+2.7 
+2.3 


+4.7 
-0.2 


+6.3 
-2.2 


a s (M T ) 


0.334 ± 0.022 


+0.7 
+0.1 


-0.7 
+ 1.3 


-1.6 

+2.2 


{ss)/(uu} 


0.8 ±0.2 


-8.7 
+ 7.7 


-9.9 
+8.9 


-10.6 

+9.6 


fK 


113 ± 2MeV 


-1.8 
+ 1.7 


-1.4 
+ 1.4 


-1.2 
+ 1.1 


Total 




+33.6 
-44.3 


+25.0 
-29.5 


+21.3 
-23.0 



TABLE I: Central results for m a {M T ) extracted from the dif- 
ferent moments, as well as ranges for the main input param- 
eters and resulting uncertainties for m s . 

cuss these here. 

The analysis proceeds in complete analogy to our pre- 
vious work [4]. In table I, we show a detailed account of 
our results. The first row displays the values of m s (M T ) 
obtained from the different moment sum rules and cen- 
tral values for all input parameters. In the following rows, 
we have listed those input parameters which dominantly 
contribute to the uncertainty on m s , the ranges for these 
parameters used in our analysis and the resulting shift 
in m s . Only those parameters have been included in the 
list which at least for one moment yield a shift of m s 
larger than 1 MeV. Finally, in the last row, we display 
the total error that results from adding the individual 
uncertainties in quadrature. 

Taking a weighted average of the strange mass values 
obtained for the different moments we then find 

m s (M T ) = 84 ±23 MeV, 
=> m s (2 GeV) = 81 ± 22 MeV , (7) 

where the uncertainty corresponds to that of the (4,0) 
moment. The dominant theoretical uncertainties in the 
result of eq. (7) originate from higher order perturba- 
tive corrections as well as the SU(3)-breaking ratio of 
the quark condensates (ss) / (qq) [34] which arises in the 
dimension-4 contribution to eq. (4). A detailed discussion 
of all input parameters can be found in rcf. [4]. 

In our previous analysis [4] , based on the ALEPH data 
[10], it was observed that m s displayed a strong depen- 
dence on the number of the moment k, decreasing with 
increasing k, and it was speculated that this behaviour 
could be due to missing contributions in the higher en- 
ergy region of the spectrum [35] . With the recent CLEO 
and OPAL data [21, 22], finding a larger branching frac- 
tion of the if~7r + 7r~ mode, the decrease of m s is now 
much reduced, although still visible. This issue needs to 
be clarified further once even better data are available. 



SIMULTANEOUS FIT OF Vus AND M s 

In principle, it is also possible to perform a simulta- 
neous fit to V us and m s from a certain set of (k, I) mo- 
ments. As soon as more precise data are available, this 
will be the ultimate approach to determine V us and m a 
from hadronic r decays. With the current uncertainties 
in the data and the question about a monotonous k- 
dependence of m s , a bias could be present in the method. 
Furthermore, the correlations between different moments 
are rather strong and also have to be included on the the- 
ory side. 

Here, we shall restrict ourselves to a simplified ap- 
proach where all correlations are neglected. For the si- 
multaneous fit of V us and m s , we employ the five Rlf 
moments (0, 0) to (4, 0) which have also been used in our 
previous analysis [4]. Performing this exercise, for the 
central values we find: 

\V US \ = 0.2196, m s (2GcV) = 76MeV. (8) 

The expected uncertainties on these results should be 
smaller than the individual errors in eqs. (6) and (7), 
but only slightly since the correlations between different 
moments are rather strong. 

The general trend of the fit result can be understood 
easily. m s from the OPAL data turned out lower than 
our global average m s (2 GeV) — 95 ± 20 MeV considered 
in section 3. Thus, also the corresponding SR T , t h is lower, 
resulting in a slight reduction of V us . Furthermore, the 
moment-dependence of m s is reduced in the fit. Never- 
theless, leaving a detailed error analysis for a forthcoming 
publication [36], at present, we consider eqs. (6) and (7) 
as our central results. 

CONCLUSIONS 

Taking advantage of the strong sensitivity of the 
flavour-breaking r sum rule on the CKM matrix element 
V us , it is possible to determine V us from hadronic r decay 
data. This requires a value of the strange quark mass as 
an input which can be obtained from other sources like 
QCD sum rules or the lattice. The result for V us thus 
obtained is 

\V US \ = 0.2208 ±0.0034, (9) 

where the error is largely dominated by the experimen- 
tal uncertainty on i? Tj s, and thus should be improvable 
with the BABAR and BELLE r data sets in the near 
future. Already now, our result is competitive with the 
standard extraction of V us from K e ^ decays [37-40] and a 
new determination from fx/ fn as extracted from the lat- 
tice [41, 42]. The resulting deviation from CKM unitarity 
then is 

1 - \V ud \ 2 - \V US \ 2 - \V ub \ 2 = (2.9 ± 1.8) • 10- 3 , (10) 
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being consistent with unitarity at the 1.6 a level. 

For the strange mass determination, we have used the 
three moments (2,0) to (4,0), with the result 

m s (2GcV) = 81±22MeV. (11) 

Our value for m s is on the low side of previous strange 
mass determinations, but certainly compatible with 
them. It is also on the borderline of being compatible 
with lower bounds on m s from sum rules [16, 43-46]. 

Finally, we have performed a simultaneous fit of V us 
and m s to the five moments (0,0) to (4,0). Our central 
values arc completely compatible with the central results 
of eqs. (6) and (7). Anticipating a detailed analysis of the 
correlations between different moments, these findings 
should be considered as an indication of the prospects 
for the future when more precise experimental data will 
become available. 
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